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In a recent work Killip and Nenciu gave random recurrences for the characteristic polyno- 
mials of certain unitary and real orthogonal upper Hessenberg matrices. The corresponding 
eigenvalue p.d.f.'s are /^-generalizations of the classical groups. Left open was the direct 
calculation of certain Jacobians. We provide the sought direct calculation. Furthermore, 
we show how a multiplicative rank 1 perturbation of the unitary Hessenberg matrices pro- 
. vides a joint eigenvalue p.d.f. generalizing the circular /3-ensemble, and we show how this 

joint density is related to known inter-relations between circular ensembles. Projecting the 
|**y^ joint density onto the real line leads to the derivation of a random three-term recurrence for 

£Lh ' polynomials with zeros distributed according to the circular Jacobi /3-ensemble 

1 Introduction 

Consider the classical group U(N) of NxN unitary matrices. There is a unique measure cLh U — the Haar 
measure — which is invariant under both left and right multiplication by a fixed unitary matrix, thus 
| giving a uniform distribution on the group. The corresponding eigenvalue probability density function 

JO ■ (p.d.f.) has the explicit form (see e.g. [5]) 

)Q ■ TT | e ^-e^l 2 fill 

l/j ■ v ' l<j<fc<iV 
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and this in turn is of fundamental importance in recent applications of random matrix theory to combi- 
natorial models [201 Bj , analytic number theory ^H] and the quantum many body problem [H] . 

A basic question is how to best sample from p.ljl . Until very recently, the only method available 
has been to first generate a member of U(N) according to the Haar measure, by for example applying 
the Gram-Schmidt orthogonalization procedure to the columns of an N x N complex Gaussian matrix, 
then to calculate the eigenvalues of the resulting matrix. However, inspired by recent work of Dumitriu 
and Edelman this situation has been dramatically improved upon by Killip and Nenciu [TJ\. Thus 
augmenting ideas from jS] with results from the theory of orthogonal polynomials on the unit circle, these 
authors have provided an explicit unitary Hessenberg matrix, with positive elements on the subdiagonal, 
which has for its eigenvalue p.d.f. the /3-generation of 

In general the characteristic polynomial xjv(A) of such matrices can be calculated from the coupled 
recurrences 

Xfe(A) = Ay fe _i(A) - a fc _iy fe _i(A) 

Xfe(A) = Xfc-iO) - Aa fc _iy fe _i(A) (1.3) 
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(k = 1, . . . , N) where Xo(A) = Xo(A) = 1 and furthermore Xfe(A) = A fc Xfc(l/A). For the unitary Hessenberg 
matrix relating to (|1.2|) . the parameters {£*}}j=o,...,iV— l are random variables with distributions specified 
in ^7] (see l|4.3|) below). As a consequence of this result the joint distribution (|l.lfl . or more generally 
l|1.2fl . can be sampled by simply iterating to generate xnW, then computing its roots. 
The problem of efficiently sampling from the p.d.f. 



N 



C N {a,b;P) 



J]|l-e i9 f Q + 1 |l + e ' 9 f b + 1 J] | e *i -e^fll-e^+^f, (0 < fl, < tt) (1.4) 



was solved according to the same strategy in |17j . Here the underlying unitary Hessenberg matrix is 
real orthogonal with determinant +1, and thus the characteristic polynomial X27v(A) has real coefficients. 
In this case Xfc(A) = A fc Xfc(l/A) and so only the first of the recurrences in is required. Note 

that the eigenvalues of a real orthogonal matrix with determinant +1 come in complex conjugate pairs 
e ±w ; (|1.4|l is the joint distribution of those with < 8 < n. The case (3 = 2, (a,b) = (±|,±i) (the 
signs chosen appropriately) of (|1.4|l gives the eigenvalue p.d.f. for matrices from the real orthogonal and 
symplectic classical groups with Haar measure (see e.g. |S]). Like their counterparts from U(N), such 
random matrices are of fundamental importance in applications of random matrix theory to combinatorial 
models [201 01 , analytic number theory jT^J and the quantum many body problem JOj • 

In the work ^7], Killip and Nenciu left open two questions concerning the direct computation of 
certain Jacobians, one relating to unitary Hessenberg matrices corresponding to (|1.2|l . and the other to 
real orthogonal Hessenberg matrices corresponding to l|1.4|l . Earlier, Dumitriu and Edelman |H] had left 
open an analogous question in the case of tridiagonal matrices corresponding to the Gaussian /3-ensemble 
p.d.f. 



i N 

G/3AT 



1=1 l<j<k<N 

It was remarked in [§] that one of the present authors (P JF) had communicated a direct derivation of the 
sought Jacobian for the change of variables from the elements of a tridiagonal matrix, to the eigenvalues 
and the first component of the eigenvectors. A primary purpose of this article is to show how a similar 
approach can be used to answer the two questions left open in |17| . We begin in Section 2 by presenting the 
calculation for the Jacobian in the case of a tridiagonal matrix. In Section 3 this calculation is extended 
to provide a direct calculation of Jacobians relating to unitary and real orthogonal Hessenberg matrices. 
Also shown is how portions of the working in |17| reliant on the theory of orthogonal polynomials on 
the unit circle, can alternatively be derived within a matrix setting. In Section 4, it is shown how a 
certain rank 1 multiplicative perturbation of unitary matrices leads to the derivation of a joint eigenvalue 
p.d.f. generalizing (|1.2|) . An integration formula associated with this p.d.f. is discussed, which in turn 
is shown to include as special cases known inter-relations between circular ensembles. Furthermore, the 
multiplicative perturbation is used to give an alternative derivation of these inter-relations. 

The Cayley transformation of the distributions obtained in Section 4, projecting the unit circle to the 
real line, are studied in Section 5. This leads to a random three-term recurrence for the projection onto 
the real line of polynomials with zeros distributed according to 

1 N 

j^^l-e-T n Ie*>-e*'l*, (1-5) 

v ' ' 1=1 l<j<k<N 

which with 2c = j3 is known as the circular Jacobi /3-ensemble |Hj. In the case a = this recurrence 
scheme is distinct from the scheme (|1.3f) . 
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2 Calculation of a Jacobian for tridiagonal matrices 



Let 



T = 



0-n b n -i 
b n -l a n-l °n-2 

bn-2 a n-2 &n-3 



(2.1) 



b 2 a 2 bi 
b\ <X\ 

be a general real symmetric tridiagonal matrix. The problem posed in [S] is to compute the Jacobian for 
the change of variables from the description of T in terms of its entries, to the description in terms of 
eigenvalues and variables relating to its eigenvectors. 

As is well known, and easy to see by direct substitution, for each eigenvalue A& and corresponding 
eigenvector Vk, once the 1st component —: q%, qk > 0, of Vk is specified, the other components are 
then fully determined by {Afc} and the elements of T. However only n — 1 of these components are 
independent due to the relation 

n 



fe=i 



which itself is a consequence of T being symmetric and thus orthogonally diagonalizable. Thus the In — 1 
variables 

a := (a„,a„_i, . . . ,a%), b := (b n -i, ■ ■ ■ , h) (2.3) 



can be put into 1-to-l correspondence with the 2n — 1 variables 

A := (Ai, . . . , A„), q := (91, . . . , q n -x) 



(2.4) 



where Ai > • • • > A„ and qi > 0. The Jacobian for the change of variables from l|2.3(l to (|2.4(l can be 
computed directly using the method of wedge products (for an introduction to the use of this technique 
in random matrix theory see 

We will first isolate results required in the course of the calculation. 

Proposition 1. Let (X)ij denote the ij entry of the matrix X. We have 



Also 



and 



((T-XI n )- 1 ) 11 =J2 



3=1 



A, - A 



n™- 1 h 2i 



1<2<j <n 



(let 



[Ai-^]^ v .., 2 n- 1 b-Ar i ]i= l ..... 2 »- l 



n (a, 

l<j<fc<n 



(2.5) 



(2.6) 



(2.7) 



Proof. The identity H2.5(l . which is well known, follows by writing the matrix entry as an inner product, 
and decomposing the vectors in this inner product as eigenvectors. The identity l|2.6|l is contained in 
It can be derived from l|2.5|l by using the fact that for a general n x n non-singular matrix 



(let X 



n-l 



dctX 



(2- 
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where X n _\ denotes the bottom right n — 1 x n — 1 submatrix of X, introducing the corresponding 
characteristic polynomials P n _i(A), P„(A), and making use of the three term recurrence 

P fe (A) = (A - a k )P k -i(X) - bl^Pk^iX), P (A) := 1. 

For the identity l|2.7[) . note that both sides are homogeneous symmetric polynomials of degree 1n(n — 1). 
Furthermore, the determinant and its first three derivatives with respect to Ai vanish at Ai = \i- As 
a consequence, it follows that the determinant must in fact be proportional to the fourth power of the 
product of differences as given in the r.h.s. The fact that the proportionality constant is unity follows by 
comparing coefficients of (A^A^ • • • A™ -1 ) 4 on both sides. □ 

The Jacobian can now be computed according to the following result. 
Theorem 1. The Jacobian for the change of variables fff.ffi to \2.$ is equal to 

1 K=i h 



nn 
i=i <H 

Proof. Rewriting 12. 5|) in the form 



(2.9) 



((Z n _AT)- 1 ) 11 =^ r ^ x -, (2.10) 



recalling the explicit form of T from 1|2.1|) , and equating successive powers of A on both sides gives 



n 

1 = Erf 

n 

a « = E q j X 3 
3=1 

n 

<-i - E^ 



3=1 

n 

* + b 2 n _ 2 b 2 n _ 1 = J2^ X 



4 
3 

3 = 1 



3=1 



* + a 1 bf...bl_ 2 bl_ 1 = Y.$ X T~ l - ( 2 - n ) 

3=1 

Here the * denotes terms involving only variables already having appeared on the l.h.s. of preceding 
equations. Thus the variables a n , a n —l-> b-n-2, ■ ■ ■ occur in a triangular structure. Upon taking 
differentials, the first of these equations implies 

n-l 

q n dq n = ~ ^ qjdqj. 

3 = 1 
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For the differentials of the remaining equations, we use this to substitute for dq n , and then take wedge 
products of both sides. On the l.h.s., the triangular structure gives 

n-l 

2 n ~ 1 Y[bf- 3 daAdb (2.12) 

where 

n n—1 

da :— ^ da,j, db :— A dbj. 
j=i j'=i 
On the r.h.s. the wedge product operation yields 

n-l 

2 n_1 <^ TT ?? dct \l X l - H) 3=i.--. 2 «-i [?'A( _1 ]j=i 2n-ildAA df (2.13) 

J |_ fc=l n — 1 k = l,...,n J 

where 



dX := A cLV,-, dq := A c?q 
i=i j'=i 



n-l 

'•V 



(a common factor 2^ has been removed from column k, k = 1, . . . , n — 1 of the determinant, as has a 
common factor q 2 from columns n — 1 + k, k = 1, . . . ,n). 
By definition the Jacobian J satisfies 

daAdb= JdX A cfcf. (2.14) 
Equating i|2.12[) and (|2.13() . and using (|2.7I) to evaluate the determinant then shows 



<?« 1I,-=1<& V n j = l V ' l<J<k<N 



Recalling QZlfy gives the form of J □ 

In [H] indirect methods are used to derive (|2.9() but with the factor l/q n not present. As noted in 
[T7) . the reasoning of [S] is most suited to working with the variables pj = q 2 , and doing so eliminates 
this apparent discrepancy. 

3 Calculation of a Jacobian for unitary and real orthogonal Hes- 
senberg matrices 

3.1 Preliminaries 

In general a unitary upper triangular Hessenberg matrix H = [Hi t j]ij=o,...,n-i with positive elements 
along the sub-diagonal is parametrized by n — 1 complex numbers ao, oci, . . . , a n -2 with |ay| = 1 and a 
further complex number a n -i with |a n -i| < 1. Setting a-i := — 1, pj := y/l — |ay| 2 (J = 0, . . . , n — 2), 
one can check that if the diagonal entries are specified as H^i — — ati—idti, and subdiagonal entries as 
Hi+i.i = pi, then the remaining non-zero entries are given by 

j-i 

Hij = -a^iOLj Y[pi, i< j- (3.1) 
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Let Aj = e lBj (j — l,...,n) denote the eigenvalues of H and let qj denote the modulus of the 
first component of the corresponding normalized eigenvectors (the {qj} thus satisfy (|2.2p ). With the 
{9j} ordered, there is an invertible 1-to-l correspondence with the parameters {<x/}j=o,....n-i- Our 
interest is to directly compute the Jacobian for the change of variables from {aj}j = o n—i to {#i}i=i,...,n, 

{<?i}i=l,...,n-l' 

In preparation for the derivation of this result, note that with denoting the top k x k block of H , 

Xfc(A) := det(A/ fe - H k ) 

satisfies (|1.3[1 (see e.g. J3]). Also of interest is a variant of the characteristic polynomial of the bottom 
k x k submatrix. In relation to this, note that with the involution aj i— > — a,ja n -\ (J = 0, . . . , n — 2), the 
bottom k x k submatrix, after reflection in the anti-diagonal, becomes equal to the top k x k submatrix 
but with aj i— ► a n _2— j 0=0,..., n — 2). Let the characteristic polynomial of the bottom kxk submatrix 
with the replacements aj i— > — oija n -\ (j = 0, . . . ,n — 2) be denoted x|(A)- We see that this polynomial 
satisfies the recurrence <|1.3[) with aj replaced by — a n _2-jCt„_i m (|1.3|) . 

XkW = A Xfc-i(A) + a«-i-fca„_ixt-i( A ) 

Xfc(A) = Xfc-xW+Afin-i-fcOn-ixLxW (3-2) 

(fc = 1, . . . , n) where Xo(A) = Xo(A) = 1 and x|(A) = A fe Xfe(l/A). These recurrences can be used to derive 
the analogue of (|23|) [T7| 

Proposition 2. i^or i/ie unitary Hessenberg matrix specified by and surrounding text, we have 

n ^-tf- ^y" (3.3, 

Proof. We follow the strategy sketched to prove (|2.(il) in the proof of Proposition^ Analogous to 12.51) 
we have 

n 2 

((H-XI n )- 1 ) 11 ^Y,J S ^x (3 ' 4) 
where ^ and Aj relate to H as specified below (|3.1[1 . Using l|2.8ll . (|3.4I) can be rewritten as 

Xn — 1 (A) | 

nr=i(A-Ao = ^a-a-' (3 - 5) 

which in turn implies 

1~L=1 \Xn-l{^i)\\otj^-ot :j a n --i 



IK 



, , rii<i<j<„ |A; - Aji 2 

From this we see <|3.3[1 follows if we can show 

n n — 2 

ni^-i(Ai)i=n( i -Ki a ) n_i ~ i - ( 3 - 6 ) 

i=l 1=0 

To establish g3J we will use (E2J). With A^ p) denoting the jth zero of Xp(A), it follows from 
that 

X^VAf) = i(l-|a B - fc -i| 2 )xLi(lAf) 

x|(Af) = (l-|a„_ fe _i| 2 )xt 1 (Af). (3.7) 
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Introducing the factorizations 

Xk-i( x ) 
we deduce from (|3.7|) that 



fe-i 
Ifc» 



.(k-i) 



fc 



(1 



(1 



-fc-1 



"n-fc-1 



J ) fe Il(l/A 



k-i 
( fc )\fe TT .-.6 



i=l 

fe-1 



3=1 



These latter two equations together imply 

k 



3=1 



2\2 + l 



1=0 



Making further use of the first equation in ((3.7|) . setting k — n, and noting \x[ n ^\ = 1 gives ((3.6|) 



□ 



I n ((3.3(1 is derived using a different strategy relating to the determinant of the Toeplitz matrix 
formed from the moments of the underlying measure. 

Also required is a determinant evaluation analogous to 12.711 . 



Proposition 3. We have 



dct 



K - K 



X~ J 



- j,k=l,...,n— 1 



[«Aj?]*=i 



-iK 3 



J=l>.- 

fe=l, 



ni 



<j<fe<n(Afc 



nn \2n-Z 
1=1 A l 



(3.8) 



Proof. By inspection the determinant is a symmetric function of Ai, . . . , A„ which is homogeneous of 
degree n. Upon multiplying columns 1 and columns n by A 2 ™ -3 we see that the determinant becomes a 
polynomial in Ai, so it must be of the form 



1 



2n-3 



p(Xi, . . . , A„) 



nr=i A 

where p is a symmetric polynomial of Ai, . . . , A n of degree 2n(n — 1). 

We see immediately that the determinant vanishes when Ai = A2. Furthermore, it is straightforward 
to verify that its derivatives (Ai gf-)^ (j = 1, 2, 3) also vanish when Ai = A2. The polynomial p must thus 
contain as a factor rii<j<fe<n(Afe — A.,) 4 . As this is of degree 2n(n — 1), it follows that the determinant 
must in fact be proportional to 13.8(1 . 

On the r.h.s. of ((3.8(1 . the coefficient of Yi?=i Af 2,1+3 [ s unity. In the determinant, let us add 
(n — 1) times the first column to the nth column. Then we see that the coefficient of A^ 2n+3 is given by 
the cofactor coming from multiplying together the (2n — 2, 1) and (2n — 4, n) elements. In the cofactor 
we add (n — 2) times the first column to the (n — l)st column. The coefficient of A7/ 2,i+7 is given by the 
cofactor coming from multiplying together the (2n — 3, 1) and (2n — 5, n — 1) elements. Proceeding in 
this manner we see that the coefficient of Y\x=i A 4 ^ ^ 2n+3 is also unity in the determinant. □ 

As remarked in the Introduction, the approach to the ensemble 1(1.4(1 given in 17 is via 2n x 2n 
(n = N) real orthogonal Hessenberg matrices with determinant +1. The elements being real implies 
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{ctj}j=o,...,2n-i are real, while the determinant equalling +1 implies a2 n —i = ~ !■ Thus there are 2n — 
1 independent real parameters ao, . . . , a^n— i- I n the corresponding eigen-decomposition, there are n 
independent eigenvalues Xj = e i6j (j = l,...,n, < 6j < it) and n — 1 independent variables qj 
(j = 1, . . . ,n— 1) where is the square of the first component of both the eigenvalues Xj and Xj. Left 
open in is the problem of a direct calculation of the corresponding Jacobian. For this the analogue 
of Proposition is required. 

Proposition 4. \l r /Jj For a In x 2n real orthogonal Hessenberg matrix of determinant +1, parametrized 
in terms of the real parameters {ai}i=o,...,2n-2, < 1, we have 

Il^-fl II |A I -A J f|A^l/A J f^2" 11 ^ (1 ai \ (3.9) 

i=l Al l<i<j<n llt=l«t 

n 2n — 2 n In — 2 

nil -A.-I 2 =2 n a-**). ni i+ ^i 2 = 2 ( 3 - 10 ) 

Proof. Denoting the Hessenberg matrix in question by H, the analogue of ll.'j.SH reads 

(fe-Atf>-')„ = 5 I>j(— + TT ^). (3.11) 

j = l J J 

Analogous to the reasoning underlying (|3.5|) . the l.h.s. is equal to X2n-i(l/A)/Ax2n(l/A)- We thus have 



X2n-l(Aj) = 1 

X 2 „(Ai) 2 



3? (i = l,...,2n) 



where Xj+ n = Xj, qj+ n = f?j- Taking the product over j = 1, . . . , 2n, making use of 1|3.6|1 . then taking the 
square root gives H3.9[) . For the results (|3.1()jl . one notes 

n n 

II I 1 - A,| 2 = X 2n(l), LI I 1 + = X2»(-l), 

3=1 J=l 

while from (JOJl Xfc+i(A)U=±i = (A - a fc A fe )xfe(A)| A= ±i. □ 

We remark that in ^7] 1|3.9[) is deduced by making use of (|3.3[) , which in turn is derived using formulas 
relating to the underlying measure. Our derivation of (|3.1U|) is the same as that in |17j . 
We must make note too of a further determinant evaluation. 

Proposition 5. We have 

det [[X{ + X-J (X{ + \j(Xi - A^)] n ^., 2? -/ 

n 

= n(A J -l/A,) I] (Afe-A.^a/Afe-l/A.O^A.-l/A^^l/A.-Afe) 2 . (3.12) 

■j = l l<j<k<n 

Proof. We see that the determinant is a symmetric rational function in Ai, . . . , A n , and is antisymmetric 
under the mapping Ai > 1/Ai for any i = 1, . . . , n. It must thus be of the form 

n 

n^-l/A^Ai,...,^) (3-13) 

3=1 

where q is symmetric and unchanged by the mapping Ai i— ► 1/Aj. Noting too that the determinant 
vanishes when Ai = A2, we see that q must contain as a factor 

[] (A fe -A J ) 2 (l/A fe -l/A J ) 2 (A,-l/A fe ) 2 (l/A J -A fe ) 2 . (3.14) 
l<j<k<n 
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The highest order term (in degree) of l|3.14|l multiplied by Y\j =1 (Xj — 1/Xj) is IIJ=i Y\i<j<k<n(^ k ^ -\?') 4 - 
On the other hand the highest order term in degree in the determinant is 



dct 



[A{-^]i=i,.... a »-i[jA{],-i...., a »-i 



According to (|2.7|l this evaluates to the same expression, so in fact q must be exactly equal to 13.14fl . □ 
3.2 The Jacobians 

Using a similar strategy to that used to derive the Jacobian (|2.9|) in the proof of Theorem^] the results of 
the previous subsection together with the method of wedge products allows the two Jacobians evaluated 
by indirect means in |17| to be derived directly. 

Theorem 2. Consider unitary Hessenberg matrices with entries specified by \3.1\l and surrounding text. 
The Jacobian for the change of variables from {<Xj}j=o,...,n— 1 to {0»}i=i ,n> {<7i}i=i,...,n-i is equal to 

ICfr-W) , (3.15, 
q n Ui=i Qi 

Consider 2n x In real orthogonal Hessenberg matrices as specified above Proposition^ The Jacobian 
for the change of variables from {otj}j=o in-i to {Oi}i—i t „, lV ,, {(fa}i=i,...,n-i equal to 

a"- 1 nf: ' 2 (i-hi 2 ) (316) 



? n nr=i * ulTo 2 ( i - «fe) i/2 (! + (-i) fe « fc ) i/2 ' 

Proof. In relation to i|3.15fl we begin by equating successive powers of A on both sides of (|3.4|l . Recalling 
the explicit form of H given by l|3.1fl and surrounding text this gives 

n 

i - £«? 

3=1 



0: n 



3=1 



* + "iPo = E 5 I A 



2 

3 

3=1 



* + a 2 PoPi = E 9 f A 



3 

3=1 



* + a„-ipgp?---/£_ a = (3.17) 



Z^i ^3 3 

3 = 1 

where the * denotes terms involving only variables already having appeared on the l.h.s. of the preceding 
equation. Thus as in the corresponding equations (|2.11|) in the tridiagonal case a triangular structure 
results. We know that otj (j = 0, . . . , n—2) has an independent real and imaginary part, while a n _i := e 1 ^, 
Xj := e l6j (j = 1, . . . , n) have unit modulus. Consequently the number of equations can be made equal 
to the number of variables by firstly using the first equation to eliminate q\ in the subsequent equations, 
then appending to the list the complex conjugate of all but the last of the remaining equations. 
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Let us take differentials of these 2n — 1 equations, then take wedge products of both sides. Because 
of the triangular structure, we obtain on the l.h.s. 

n-2 

P 2 oP 2 i ■ ■ ■ pI-i II Pn-l-2dct A d<f>, (3.18) 



1=1 



while this operation on the r.h.s. yields 



n-l 

i=i 



det 



k n J j,k=l,...,n-l 

[Afc _ A"]fc=l,. ..,„_! 



jA{ 



-j 



[nA^]fe=i,... ) „ 



3 = 1,.. 

fe=i, 



dO Adq 



««ru 3 n iA fc -A,i 4 ^Adq 

j=l l<j<fc<n 



(3.19) 



where the equality follows upon using the determinant evaluation (|3.8|l . 
Analogous to (|2.14|) . by definition the Jacobian J satisfies 

da Ad(f)= Jd9 A dq. 

Equating l|3.18|l and (|3.19|) and making use of l|3.3|l gives (|3.15|) . 

Consider next the derivation of (|3.16|) . Proceeding as in the derivation of (|3.17|) . expanding (|3.11|) in 
powers of A, we obtain 

n 

i - Erf 

i=i 
1 ™ 



Of) 



5 i; 



* + «2n-2Po ' ' 'Pins 



^E^^ + Af- 1 ). 



The l.h.s. again exhibits a triangular structure, and furthermore all quantities on the l.h.s. are real. 
Taking the differentials of both sides, and forming the wedge product of the l.h.s. 's of all but the first 
equation gives 



2n-3 



n P f 2n - 2 - i) da. 



(3.20) 



On the r.h.s., after substituting for q n dq n using the differential of the first equation, this same procedure 
gives 

n-l 

2 ~ n il II a? | det N + V - (K + x k j )]i=Y- 2 :-i M x * - V^t-r 2 n A I Sh d<r - (3 - 2l) 

i=i 

Here the Jacobian J satisfies 

da = Jd9 A dq, 

so <|3.2U[1 and (|3.21|) (with the determinant evaluated according to l|3.12ll ) together give a formula for J in 
terms of {eft}, {c^i} and {A,;}. The latter set of variables can be eliminated by making use of Proposition 
H and (|3~T?)|) results. □ 



10 



4 A multiplicative rank 1 perturbation of unitary matrices 



4.1 Circular analogue of the Dixon- Anderson density 

Let e\ denote the n x 1 unit vector (1,0,..., 0) T . Let t be a complex number with \t\ — 1. Then the 
matrix 

/„ - (1 - t)eief 

is a unitary matrix differing from the identity only in the top left entry which is t. Our interest in this 
section is the eigenvalue distribution of 

U:={I n -{l-t)Si^)U, (4.1) 

for U a given unitary matrix. Such multiplicative rank 1 perturbations are discussed for example in 
PP. The term multiplicative perturbation is used because U is obtained from U by multiplication of the 
first row by the unimodular complex number t, while the term rank 1 is used because the multiplicative 
perturbative factor differs from the identity by a rank 1 matrix. We will see that for U a member of the 
circular /3-ensemble, a joint eigenvalue p.d.f. generalizing i|1.2fl results. 

First a rational function having as its zeros the eigenvalues of U will be specified. 

Proposition 6. Let U be an n x n unitary matrix with distinct eigenvalues e *, . . . , e n , and denote the 
corresponding matrix of eigenvectors by V — [Vjk]j,k=i n- The eigenvalues of U as specified by \4-l\j 
occur at the zeros of the rational function 

e ie i\ih -I 2 

^) = i+(*-i)£y=l (42) 

Proof. Noting from l|4.1|l that U = U — (1 — t)U', where U' is the n x n matrix in which the first row 
is equal to the first row of U, and all other rows have all entries zero, we see 

V~ 1 UV = diag[e l91 , . . . , e 4e "] + (t - l)[ Vlj v lk e l %, k=1 ,..., n . 

Thus U has the same spectrum as a matrix which consists of a rank 1 multiplicative perturbation of a 
diagonalized unitary matrix. The characteristic polynomial of this matrix can be factorized as 



l[(e tei - A) det k- k + (t- l)v ljVlk e Wk /(e l9 > - A) 



1=1 

and the zeros must occur at the zeros of the determinant. Noting the simple determinant evaluation 



n n 1 

det [«,•$,-,*. + l]j,fc=i,..., n = TT^'f 1 + ^ ~ 



i=l ] = 1 ■> 

the sought result follows. □ 
We remark that Proposition [S] can be extended to the case that each eigenvalue e lSj has multiplicity 



Thus with v[ s ) denoting the first component of the sth independent eigenvector corresponding to 



e ie i, we replace \ Vlj \ 2 in g^J by J2"=i \ v ij\ 2 - 

Let us suppose now that the matrix U is a unitary upper triangular Hessenberg matrix parametrized 
as specified by l|3.1f) an d surrounding text. One of the main results of |17| is that the parameters 
{ctj}j=o,... t n-i can be chosen from particular probability distributions so that the eigenvalue p.d.f. of U 



11 



is given by i|1.2|) ■ The probability distributions in question are parametrized by a real number v > 1 and 
denoted by 8„. For 2/ > 1, the support of Q v is the open unit disk \z\ < 1 in the complex plane, and the 
distribution is specified by the p.d.f. 

^__i (1 _ N2) (-3)/2. 

For v = 1, the support is the unit circle \z\ — 1, and 0i denotes the uniform distribution. Proposition 
4.2 of [T7| tells us that if 

a n -j-i ~ Q/3j+i (i = 0, . . • , n - 1) (4.3) 

then the corresponding eigenvalue p.d.f. is given by (|1.2|) . Furthermore, it tells us that the modulus 
squared of the first component of the eigenvectors |i>ij| :— fij have the distribution with measure 



T^rK 72 - 1 ^, (4.4) 

C/3N X , 



where 

r"(/?/2) 

i^JV — Y(f3N/2) ' ' = ' ' fJ " n ~ 1 ' 

This is an example of the Dirichlet distribution. 

The latter fact motivates the study of the zeros of l|4.2|) with the \vij\ 2 distributed according to 
the Dirichlet distribution. We will find that a conditional p.d.f. relating to i|1.2[l results provided the 
distribution of t is appropriately chosen. First, some preliminary results must be established. 

Lemma 1. Suppose in \4-^V ^at 

< 6>i < 2 < ■ ■ ■ < n < 2tt. (4.5) 
The function C(A) has exactly n zeros occurring atX= e^ 1 ,. . .,e^», where 

i _l<Vi<fi (< = l,...,n, 6» :-6» n mod27r). (4.6) 
Furthermore, with Xj :— e ' , \j := e 1 ^ , we /iawe 

-(t-l)A^ = n^r^-AO (i = 1) ... )n) (4J) 

n n 

i=i i=i 

Proof. The fact that there are exactly n zeros of unit modulus follows from the relationship of C n (A) 
to the characteristic polynomial of a unitary matrix. The interlacing property is well known pp. It can 
be seen graphically by writing (|4.2|l in the form 



C n (X) = i-^-i (cot | - H COt ( 



i=l 

where we have set t :— e 1 ^, A :— e 1 ^. 

With the zeros so identified, and the poles as evident from 14.21) . it follows that 

where use has also been made of the property C n (A) — * 1 as A — > oo. Comparing residues in l|4.2|l and 
H4.9I I gives (I4.7|l . while setting A = gives (|4.8|l . □ 

The Jacobians for some change of variables are also required. 
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Lemma 2. Let J be the Jacobian for the change of variables {qj}j=i,..., n -i U {</>} to {ipj}j=i,...,n- We 
have 



J==|l-t|-("- 1 ) [J 

l<j<fc<n 

Proof. By definition J is positive and satisfies 

dq f\dt — Jdtp. 

Now 

dqAd(/> = (t - f - f )q A d(/) 



Afe — A, 



(4.10) 



(4.11) 



det 




\d(t-i)qr 




r at - 








- d\, - 


J=l,...,n 
I = 1 , . . . , n — 1 




j=l,...,n_ 



d\n 



But according to l|4.7|l and l|4.8|l 



d(t - l)qi (t - St _ t 



9Aj A/ — Aj 

and so with A„ = (temporarily as a notational convenience) 

n— l 



dXj Xj 



dqh # = t TT (— gj) det 

i=i 



1 



Aj - A; 



dA. 



j,l=l,...,n 



(4.12) 



Since J is positive and satisfies (|4.11l) . it must be equal to the modulus of the terms multiplying dX in 
this expression. Evaluating the determinant as a Cauchy double alternant, and evaluating nl^Ti 1l using 
H4.7fl gives the stated result. □ 

The results of the above two lemmas allow a change of variables to be made from {qj} j=i,...,n— l U {t} 
to {A{}{=i ) ..., n . 

Theorem 3. With \vij\ 2 — qj (j = 1, ...,n) in let {qj} have the Dirichlet distribution with 



r((n - l)d + do) 



(r(d))«-ir(*) 

Further, let the parameter t in \4-%ty be determined by the p.d.f. with measure 

r 2 (i(d + (n-i)d + i)) _ do+(n _ 1)d _ 1 

2^r((n-l)d + d ) 1 1 ^ 
XTie conditional p.d.f. of {Xj — e l ^ j }j = i.... yn , given {Xj = e i 3 }j=i,...,n, is egwed to 

nr=i - i^- 1 n;=i nr=i - ^ I"- 1 



riL"! 1 - e^o+d-i Ui< j<k <n-i \ el9k ~ ^ 2<l ' 



t n 



(4.13) 



(4.14) 



(4.15) 



l<j'<fc<n 



A :-- 



r((n - l)d + d ) T 2 (i(d + (n - l)d + 1)) 



(r(d))"" 1 r(d ) 27rr((n-l)d + d ) ' 

Proof. Our proof, which at a technical level proceeds by making use of the results of Lemmas ^ and 
[21 is based on a strategy adopted for an analogous problem with real roots by Anderson |2j, and many 
years before by Dixon [H] . 
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We must change variables in the product of (|4.13|l . i|4.14|l . and the Jacobian J. We know the latter 
is given by (|4.1U|) . To change variables in (|4.13|) we use (|4.8|) . which gives 

Multiplying l|4.14ll . (|4.1UI) and (|4.1t)|) we see that the dependence on t cancels, and the expression (|4.15l) 
results. □ 



4.2 Properties of the corresponding joint density 

Let us suppose {Oj}j=i,..., n , assumed ordered as in 14. 6|) and with 6 n fixed, are distributed according to 
the p.d.f. 



I — 1 II n ~ 1 

y 1 tt TT \JQn_J6i lai+do + d-l 



TT | e »e„_ e iei|oi+rfo+d-i J~| 



(27r)™- 1 M n _ 1 ((ai+d + d-l)/2,(ai+d + rf-l)/2,d) |= , , . ; , 

M N (a,b,\) =11 . , (4.17) 

7=o ( Aj + a + 1 ) r ( Aj + & + 1 ) r ( 1 + A ) 

(for a discussion of this p.d.f. see Multiplying this with (|4.15|) gives the joint p.d.f. 

n (n,n-l)( j n\ . = A ( n ~ !) ! TT \ A6 n _ iifina 

W,> ' (2n) n - 1 M n ^ 1 ({a + a 1 +d)/2,(a + a 1 +d)/2,d)f- = Y 1 

n— 1 n — 1 n 

x J] |e^ fc - | Y[ \e l9n ~ e ie > \ ai J] \e l9k - e l9 > | J] JJ |e^ - e #1 | d "^4.18) 

l<j"<fe<n i=l l<j<fc<n— 1 3=1 J=l 

where do — 1 =: a. 

The case of (I4.18|) relevant to the circular /3-ensemble of Killip and Nenciu is a = d — 1, a\ = 1 and 
d = /3/2. Then Ij4.18|l is symmetric in {0i}j=i,...,n and in {ipi}i=i,.„, n , and may again be considered as 
variable f l|4.18|) should then be multiplied by n/27r to get the correct normalization). It corresponds to the 
joint eigenvalue p.d.f. of a unitary Hessenberg matrix with parameters distributed according to (|4.3|l . and 
thus with eigenvalue p.d.f. I|1.2|) . and the same unitary Hessenberg matrix perturbed by multiplication of 
the first row by t. The factor t is to be distributed according to (|4 . 1 4|) with do = d = (3/2. We know that 
the p.d.f. for {ft}i=i,..., n can be sampled by computing the zeros of x«(A) as calculated from l|1.4[> with 
{ctj}j=o,...,n-i chosen as specified by Q4.8|l . To sample from {ipi}i = i,..., n -i in the joint p.d.f., with the 
same {oij}j=o,... y n-i we again compute Xn(^) from the recurrences 1)1. 4f> . but now with Xo(^) — Xo = t. 

Next, let us turn our attention to integration formulas associated with (|4.18() . Since (|4.15|) is a 
conditional p.d.f. we must have 

R (2-K) n 1 Af„_i((a + ai + d)/2, (a + ai + dj/2, d) 

n—1 

x Y[ \e lB ~ - e lBl \ a + a i+ d Yl W 6k ~ \ 2d (4.20) 

1=1 \<j<k<n-l 

where R denotes the region specified by the inequalities 14.5|l . Special cases of (14.19(1 are two classical 
inter-relations between circular ensembles |7lll8| (for an extensive study of such formulas in random matrix 
theory see [TT], and for application of the Dixon- Anderson density to the cases with real eigenvalues see 

E2D- 
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Proposition 7. Let COE n , CUE n , CSE n — the circular ensembles with orthogonal, unitary and sym- 
plectic symmetry respectively — refer to the eigenvalue p.d.f. with j3 = 1,2,4 respectively. Let alt 

refer to the operation of integrating over every second eigenvalue. Let COE„ U COE n denote the p. d.f. of 
2n eigenvalues which results from superimposing two independent sequences of n eigenvalues each with a 
COE n distribution. One has 

alt(COE„ U COE„) = CUE„ (4.21) 
alt(COE 2 „) = CSE„ (4.22) 

Proof. For the first identity we require the fact that ^3] 

COE„ U COE„ oc J] \e l02k - e'^He* 8 *- 1 - e ie ^~ l \. (4.23) 

l<j<k<n 

We then see that Ij4.19|) with a% = d = 1 is equivalent to the first identity. The second identity is 
immediately seen to correspond to (|4.19|l with a-y = 2, d = 2. □ 

4.3 Matrix theoretic derivation of the COE, CUE, CSE inter-relations 

The inter-relations 14.21J1 , (|4.22(l were originally proved by establishing the same integration formulas as 
those noted in the proof of Proposition In this subsection it will be shown how random matrices can 
be constructed in such a way that both 14.21fl and (|4.22l) are immediate. 

Let us consider first 1)4. 21|) . This requires a different random matrix realization of the joint p.d.f. (|4.18(l 
in the case d = oi = l,o=0to that given in the paragraph below l|4.18jl . The theory underlying the 
construction is the following. 

Proposition 8. Let Mi be a 2n x 2n unitary matrix with real elements, which has a doubly degenerate 
spectrum with the independent eigenvalues distributed as CUE n . Let the matrix of eigenvectors be V = 
["ij'liii=i,-,2n) an d suppose the joint distribution of [ij := (i>i2j— l) + (^I2j) 2 (j — 1, is equal to 

the Dirichlet distribution {4-4^ with (3 — 2. Form the matrix M[ by multiplying the first row of Mi by the 
complex number t, \t\ — 1, where t has distribution \4-14\l with do = d = 1. Then the perturbed matrix 
M[ has for its eigenvalue p.d.f. \4- 18\) with a = 0,ai = d= l. 

Proof. Let {e iei }i = i^... in denote the independent eigenvalues of Mi. Proceeding as in the derivation of 
l|4.2J) shows that the characteristic polynomial of the perturbed matrix is equal to 

n n ,ig 

n(^-A) 2 (i+(*-DE^r)- ^ 

J'=l J'=l 

Thus M[ has n eigenvalues at {e l }i=i,..., n , an d n eigenvalues given by the zeros of the rational function 
factor in (|4.24() . We are given that the former have p.d.f. CUE„, while Theorem [3] tells us that the 
latter have conditional p.d.f. (|4.15|l with do = d = 1. Multiplying these together gives the stated joint 
distribution. □ 

To realize the matrix Mi , we begin with an element of U (n) chosen according to the Haar measure. 
This gives the eigenvalue p.d.f. CUE„, with the eigenvectors such that the fj,j := \vij\ 2 (j = 1, . . . ,n) 
have joint distribution l|4.4|l with = 2. To obtain a doubly degenerate spectrum, each element x + iy is 
replaced by its 2 x 2 real matrix representation 

x y 

-y x 
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so forming a 2n x 2n matrix with real entries. Since the corresponding perturbed matrix M[ retains all 
distinct eigenvalues of M\ , 

unpert(Af{) = CUE„, (4.25) 

where the l.h.s. denotes the eigenvalue p.d.f. of M[ integrated over the perturbed eigenvalues. On the 
other hand Proposition [3] together with (|4.23() tell us that with reference to the eigenvalue p.d.f., M[ = 
COE„ U COE„. Thus we have a matrix theoretic understanding of (|4.21l) in the sense that its validity is 
a consequence of spectral properties of M[ which avoid the need for explicit integration of the eigenvalue 
p.d.f. 

We seek a similar understanding of (|4.22() . For this we must identify an ensemble of random matrices 
with a doubly degenerate spectrum, and their perturbations, which give rise to (|4.18() in the case d = 2, 
a = a± = 1. In fact the very definition of the circular symplectic ensemble involves matrices with a doubly 
degenerate spectrum (see e.g. |H]). Thus, if for any 2n x 2n matrix X we set 

X D := Z 2n X T Z 2n , where Z 2n := I n <g> 

and select U € U(2n) with Haar measure, then matrices of the form U D U make up the circular symplectic 
ensemble. Such matrices have a doubly degenerate spectrum, and the n independent eigenvalues are 
distributed according to CSE„. Furthermore, with the matrix of eigenvectors denoted V — [vij\ij = i t ... t2ni 
one has that the fij := |i>i2j-i| 2 + |wi2j| 2 (j = l,...,n) are distributed according to the Dirichlct 
distribution (|4.4|l with f3 = 4. Consideration of these facts, together with reasoning analogous to that 
used in the proof of Proposition |H1 gives the sought realization. 

Proposition 9. Let M 2 be a member of the circular symplectic ensemble as specified above. Form the 
matrix M 2 by multiplying the first row of M 2 by the complex number t, \t\ — 1, where t has distribution 
j4-14{ l with do = d = 2. The joint eigenvalue p.d.f. of M 2 is then given by \4. ltty with a = a\ = 1, d = 2. 

Analogous to (|4.25|) it is immediate that 

unpert(M^) = CSE„. 

Because Proposition El tells us that M 2 has a joint distribution formally equivalent to COE2„, (|4.22|) is 
reclaimed as a matrix theoretic identity. 
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5 Cayley transformation 

5.1 Cauchy analogue of the Dixon- Anderson density 

In general a unitary matrix U is transformed to an Hermitian matrix H by the Cayley transformation 

In - U 

H = l Y N Tu- (fU) 

At the level of the eigenvalues, the change of variables l|5.1|l in the workings of Sections 14. II and 14.21 leads 
to a joint p.d.f. on interlacing variables on the real line, relating to the so called (generalized) Cauchy 
ensemble |21l l4"]. Properties of this allow a random three term recurrence to be derived for the (projected) 
characteristic polynomial associated with the p.d.f. (|1.5|) . 

First we apply the change of variables implied by l|5.1|l to 14.211 with the l.h.s. written as (|4.9|l . 
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Proposition 10. Consider the rational function \4-l$i) with the lower terminal of summation extended to 
j = 0. Substitute for the l.h.s. I4.i\ ) with the lower terminals in the products extended to j = 0. Under 
the change of variables 

-\) — r, Xj = ^i (j^O) 
J Xj+i 3 yj+i w 1 

z — i c — i . . 

A= , t= (5.2) 

z + t c+t v ; 

and with X — 1 we obtain 

]Xj=o( z ~ x j) = z-c (gj/go) 

+ 1) n?=i(* - i/i) 9o(^ 2 + 1) z - Vj ' 1 • j 

where 

x > Vi > xi > ■ ■ ■ > y„ > x n+ i. (5.4) 



Proof. This follows from direct substitution, together with the formula 

qo _ IliLiCfJ + i ) 



(5.5) 



which is a consequence of (I4.7|) . with lower product terminals extended to I = 0, in the case j = 0. □ 

Theorem 4. Consider the rational function $5.3p . Let {qj}j=Q,..., n -i have the Dirichlet distribution with 
measure 

r(E"= d i) A 



Let c /icsue i/ie generalized Cauchy distribution with measure 



ibutio'i 

r( 7 )r( 7 ) 



(l + ic)- 7 (l-icr 7 cfc (5.7) 



7r2 2(i-Rc 7 )r(2Re7- 1) 
w/iere 

n 

^d 4 + l = 2Re7. (5.8) 
We /iaue that the conditional p. d.f o/{^j}j=o,...,n given {yj}j=i,...,n * s e<?waZ to 

n n 

ill 1 ^) _7 ( 1 - II' 1 + (! " ^O^' 

n n 

xiin^--^- 1 n i^-^i 1 ^^ n ^ 

j=li=0 l<j</s<n 0<j<k<n 

where 

A= mm I f5 io) 

^a-Rc,) r(2Rc7 _ ! _ E » =i dj) n » =i r(d . } ■ 

Proof. The task is to change variables in the wedge product of (|5.6[) and l|5.7|) to {xj}j=o,...,n- We have 
(l+ic)-^{l-ic)-^ = |(l-zc)-^| 2 



-2R07 

% 



n;Li(i-%) 
n^=o( i_ia: 



-2Rc 7 nLiC 1 + »y0 7 ( i - %r /, „>, 

90 nr= (i+^)-(i-^ ( j 
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where the second equality follows from l|5.5|) , and the final equality uses the fact that since xi , y\ interlace 
according to (|5.4|) , 



log 

Also, for j = 1, . . . , n 



n^a+iyj) i = 2^ lo g( 1 +^)-z^ 1 °g( 1 +%)- 



J=0 3=1 



J + 1^=1,^1%-^ I 

It remains to change variables in dq A dc. Since is related to Xj and c to t — e 1 ^ as given in l|5.2[l . 

J " 2 
dq t\dc = — 1 1 + ic\ 2 TT - ; — — ; — rdx 

where J is the Jacobian l|4.10|l (appropriately modified to account for the lower terminal being 0). In 
terms of the change of variables l|5.2(l the latter reads 

j _ „ rio<j<fc<n \ X 3 ~ X k\ 

ni< J<fe < n \vj-vk\ 

and thus we have 

t -* * i n+2 1 rio<j<fe<n Fi - ^fcl /k 1q1 

dqAdc = q ^- ; — pr - J — = — j r dx. (5.13) 

|l + *^| 2 Ili^<fc<nltfi-»fcl 

Multiplying l|5.11l) . I|5.12|l and l|5.13[) gives the stated result. □ 

We remark that the conditional p.d.f. (|5.9|l appears in ^] as a generalization of a conditional p.d.f. due 
to Dixon and Anderson (see l|5.26|l below). We remark too that the distribution l|5.7(l in the case 7 real 
is the classical i-distribution. 

Integrating (|5.9|l over {xj}j=o,...,n within the region l|5.4|l we must get unity. Using this allows us to 
derive for the multi-dimensional integral 

/>OG />00 n 

I n ( T) d) = - dxi-- dx n Y[{l + ixi)-"<(l-ixi)-^ J] \ x J~ x k\ 2d (5.14) 

" / " ° 1=1 l<j<k<n 

a recurrence analogous to that obtained by Anderson for the Selberg integral. Moreover, the interme- 
diate workings will allow us to deduce a random three term recurrence for the characteristic polynomial 
associated with the p.d.f. (|1.5fl . 

Corollary 1. We have 

r ^ .2-2Rc 7 r(2Re 7 - nd - l)r((n + 

1,1+1(7, d) = tt2 r(d)|r(7)| 2 U\rf-a\d). (5.15) 

Proof. Let us denote the region 15. 4|) by R' . As remarked, integrating 1)5. 9fl over {xj},, = o,...,n within 
R' must give unity. Setting 

d\ = ■ ■ ■ = d n = d (5.16) 

this implies 



r n n n 

A d / dx ■■■dx n + iXj)-~<(l - ix^ \[ 11 \y.j - xif- 1 \ Xj - . n .-\ 

7=0 i=ll=0 0<i<k 



j=l 1=0 0<j<k<n 



n(i+^r 7+d (i-^)^ +ti n i%-^i 2d_i (s- 17 ) 

j—l l<j <k<n 
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where A d := A\ dl= ... dn=d . Thus 

If - 
-=-I n (j - d;d) = / dx ---dx n dyi---dy n T\(l+ix j y>(l-ix j y 

A d J Ri ._„ 



3=0 
n n 

xnni^-^i"" 1 n i^-wi n i^-^i- ( 5 - ig ) 

j=lj=0 l<j<fc<n 0<j'<fe<n 



On the other hand, we know from E] that 



/n n 
dyt--dy n J| Ij/j-y^inil 
l<j<fe<n j=li=0 

(r(d))" +1 



— n (5.19) 



u ' ' 0<j<k<n 



so the r.h.s. of (|5.18() is also equal to 



I n +i{r,d). (5.20) 



T{(n + l)dY 

Equating {S~T%)> and (|5~SUf) gives IpTIS) . □ 
Iterating l|5.15|l with I Q (j,d) = 1 reclaims the gamma function evaluation [8] 

n!7„(7; d) = 2 dn ("- 1 )- 2 ( Rc ^- 1 )7r"M n ( 7 - d(n - 1) - 1,7- d(n - 1) - 1, d) (5.21) 

where M N (a, b, A) is given by (|4.17|) . 

5.2 A random three term recurrence 

Consider the rational function (|5.3|) . Suppose {g^jj^o « have the Dirichlet distribution l|5.6|) with equal 

parameters (I5.16|) . and suppose c has the distribution l|5.7[) . Suppose furthermore that {yj}j=i,...,n have 
distribution with measure 

1 - 

' ' j = l l<j<fc<n 

The marginal distribution of {xj} is then given by multiplying this with (|5.9J) and integrating {yj} over 
the region R 1 \5A\ . Using Ij5.19|l gives 

1 " 

— — — l[(l + l x 3 )-^l~tx J )^ [] \^~^\ 2d . (5.23) 

Hence with p„ + i(z; 7; g?) denoting the random monic polynomial of degree n+ 1 with zeros at {aij}_y=o,...,n 
having distribution l|5.23[) . we see that (|5.3() can be written 

p n+ i(z; T ,d) _ z-c _ fa /go) , g 24 - 

(z 2 + l)p„(z; 7 -d;rf) g (z 2 + l) £j « - iy " 

A companion identity to l|5.24|l is also required. For this purpose we introduce the random rational 
function 

nn—lr \ n 

^f-^ -E^- (^5) 
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where {fJ-j} have Dirichlet distribution 

T(nd) 



We know from the work of Dixon |S] and Anderson [2] that the conditional p.d.f. of {uk} given {yj} is 
equal to 

T(nd) Ui< j<k <n-i( u j ~ u k) 



(r(d)) n Ui< j<k < n (yj-yk) 



2d- 



i=l k=l 



provided 



J/l > til > • • • > Z/n-1 > > 2/n- (5-27) 

It follows from (|5.26|l that if {yj} have distribution l|5.22[l . then the marginal distribution of {uj} is 



equal to 



T(nd) 



Yl (uj - Uk) dyi ■■ -dy n 



(T(d)) n I n h-d: d) 

V K " n W ' ; l<j<k<n-l 



3=1 3=1 k=l 



where R is the region (|5.27|) . According to (|5.17|) this can be evaluated as 

n — 1 

We therefore conclude that l|5.25[l can be written 

Pn-i(z;-f - 2d; d) = ^ ^ 
p n (z;j - d;d) ~{ z -yj' 

Comparison of l|5.24|) and l|5.28[) implies {p n {z; 7+ (n— l)d; d)} satisfy a random three term recurrence. 

Theorem 5. With B[a,f3] denoting the classical beta distribution, let 

b n ~ B[2Rej + nd~ l,nd] (n^O), b = l, (5.29) 

and let c n have the Cauchy distribution 

r( 7 + nd)T(j + nd) ic \-^+ nd Hl ~ icH^ f5 30) 



(this is \5.T\) with 7 t— > 7 + nd,). We /iave i/iai for n = 0, 1, 

] - Pn (z; 7 + (?i - l)d; d) + (l - (1 + z 2 )Pn-i(z; 7 + (n - 2)d; d), (5.31) 



Pn+i(^;7 + n d; d) 



where po := 1. 

Proof. In H5.24fl and (|5.28|l replace 7 1— > 7 + nd. We know that in general if {dj}j—o,..., n have Dirichlet 
distribution (|5.6|) . then each d, has beta distribution B[dj, X^_ ; / . dj]. Using this fact it follows that in 
<|5.24|) we now have 

qj ~ B[d, (n - l)d] (j 7^ 0), g ~ B[2Rc 7 + - 1, nd], 
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where in deriving the former use has also been made of l|5.8[l . while in (|5.28(l 



fij ~ B[d, {n-l)d\. 

The quantities are constrained by ™_ qj — 1, X}j=i A*j = 1- Substituting (|5.31|) in (|5.24|) we thus see 
that (|5.28[1 results, thereby verifying the correctness of (|5.31|) . □ 
To relate this to the circular Jacobi /3-ensemble Q1.50 . we note that with 

1 - e i9 i 

x i = h I W- ^ °)> x o = 0, 
1 + e 3 

the p.d.f. I|5.23[l with 7 1— » 7 + 2d (7 real) extended to a measure via the multiplication by dx\ ■ ■ ■ dx n , 
becomes equal to <|1.5fl with a — 2j — 2 and extended to a measure via the multiplication by d0\ ■ ■ ■ d9 n . 
Thus the zeros of the polynomial p n (z; 7 + (n — l)d; d), with 7 real, Xi, . . . ,x n say, under the mapping 

X -!Zl,=^ (j = l,...,n) (5.32) 

Xj + 1 

give for {9j} the distribution i|1.5fl with a — 2j — 2. 

As an illustration, let us consider the case 7 = 1, d = 1, which relates to averaging over U(N). There 
are a number of averages over U(N) which are known analytically. For example, with p £ Z>q, 



\TrUP\ 2 ) U( ,u(N) 



p, < p < N 
N, p>N 



2 



U(N) 



Since 

(|TrC72f) aet r W = (|5>* 

we can compute the Monte Carlo approximation 

1 M N (k) _ . 

(iTr^^^-ElE^) I (5-33) 

fc=l j=l x j T « V i« 

where use has been made of ()5.32(l and Xj refers to the jth generation of pn{z; N; 1) from l|5.31(l . 
For the latter task, we read off from l|5.29|l that 

b n ~ B[l + n, n]. 

Also, by definition the Student ^-distribution T v say, has p.d.f. proportional to (1 + t 2 /v)^^ u+1 ^ 2 so 



—T 



u=2n+l 



Significantly, the zeros of the lower order polynomials in the sequence {pj(z;j; l)}j=o,i,...,iV themselves 
allow us, via (|5.32(l . to sample from U(j). Hence l|5.33|l can be calculated for all values of N less than the 
sought value within the same calculation. Monte Carlo results obtained this way are presented in Table 
15.21 The consistency of these results is evident. 
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p\N 


2 


3 


4 
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0.98 


0.99 


0.99 


0.98 
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